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ABSTRACT 

The discovery of extended, approximately spherical weak shock waves in the hot intercluster gas in 
Perseus and Virgo has precipitated the notion that these waves may be the primary heating process 
that explains why so little gas cools to low temperatures. This type of heating has received additional 
support from recent gasdynamical models. We show here that outward propagating, dissipating 
waves deposit most of their energy near the center of the cluster atmosphere. Consequently, if the 
gas is heated by (intermittent) weak shocks for several Gyrs, the gas within 30-50 kpc is heated to 
temperatures that far exceed observed values. This heating can be avoided if dissipating shocks are 
sufficiently infrequent or weak so as not to be the primary source of global heating. Local PV and 
viscous heating associated with newly formed X-ray cavities are likely to be small, which is consistent 
with the low gas temperatures generally observed near the centers of groups and clusters where the 
cavities are located. 

Subject headings: cooling flows — galaxies: elliptical and lenticular, CD — galaxies: active - X-rays: 
galaxies - galaxies: clusters: general - X-rays: galaxies: clusters 



1. INTRODUCTION 

Successful gasdynamical models of the hot virialized 
gas in galaxy groups and clusters must satisfy several 
cardinal requirements: 

CR1 The rate that gas cools to low tem- 
peratures must be less than ~ 10 percent of 
that predicted by conventional cooling flows, 
M ~ L x / (5kT/2fj,m p ), where L x is the bolo- 
metric X-ray luminosity and T is the mean 
gas temperature. 

CR2 The gas temperature profile T(r) must 
increase from the center out to ^0.1-0.3 of 
the virial radius, i.e. T(r) must look like a 
conventional cooling flow. 

CR3 The radial abundance profiles in the hot 
gas must agree with observation and in par- 
ticular a mass 10 s — 10 9 M & of iron should be 
concentrated within ^100 kpc of the central 
E or cD galaxy. 

CR4 It is necessary to identify the source of 
heating that reduces the baryon mass fraction 
below the cosmic value in clusters/groups 
with mean temperatures less than about 3 
keV. 

The first of these requirements is necessitated by the 
well-documented absence of cooling gas in the X-ray 
spectra of galaxy groups and clusters (e.g. Peterson et 
al. 2001; Bohringer et al. 2002). Evidently some addi- 
tional source of heating is required. In view of the limited 
energy available in supernova explosions and the uncer- 
tain efficiency of thermal conduction (Voigt & Fabian 
2004), it is generally thought that the cluster gas is 



heated directly or indirectly by jets or energetic out- 
flows originating near supermassive black holes (AGN) 
in the cores of cluster-centered elliptical galaxies. But 
there is no general consensus regarding the physical pro- 
cesses responsible for the heating or the means by which 
energy is delivered to gas at large distances from the 
central AGN. In any case, the X-ray emitting gas clos- 
est to the central AGN must not be strongly heated 
since that would violate CR2. Positive temperature 
gradients within 0.1 — 0.3r„i r are typical among cool- 
ing flow clusters and groups (Allen, Schmidt & Fabian 
2001; De Grandi & Molendi 2002), although the temper- 
ature profile in a few clusters is more nearly isothermal 
(e.g. AWM-4; O'Sullivan et al. 2005). Assuming that 
the hot cluster gas is in approximate hydrostatic equilib- 
rium in an NFW mass distribution, gasdynamical models 
that produce correct temperature profiles will also have 
correct density and entropy profiles. 

Except possibly for the very center, the gas phase iron 
abundance in groups and clusters generally decreases 
from ~solar to ~ 0.3 ± 0.1 solar at the outer limit of 
observations. In cool-core clusters a significant fraction 
of the total gas phase iron mass, Mp e « 10 8 — 10 9 AT Q , 
is concentrated within ~ 100 kpc of the central E or cD 
galaxy (de Grandi et al. 2004). Metal abundance gra- 
dients in the hot gas retain an integrated record of its 
past association with stars and supernovae that must be 
reproduced by any successful model (CR3). 

To explore the gasdynamical consequences of heating, 
Brighenti & Mathews (2002a; 2003) considered a wide 
variety of ad hoc heating scenarios and described the ef- 
fect of each on the hot gas over many Gyrs. In these 
models the gas was heated cither continuously or inter- 
mittently, either symmetrically near the center or at lo- 
cations off center. X-ray cavities formed in the heated 
gas generated weak shock waves similar to those found 
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in Perseus. We concluded at that time that none of the 
many types of heating scenarios satisfied all cardinal re- 
quirements listed above. For example, if the heating is 
symmetric around the center of the flow at a level suf- 
ficient to quench the radiative cooling there, we found 
that cooling still occurred at larger radii in the flow (vi- 
olating CR1) or the radial temperature gradients in the 
heated region were invariably negative (violating CR2). 
To fully shut down a cooling flow it is necessary to heat 
the gas nearly to the ( "cooling" ) radius r coo i , the radius 
at which the gas would cool in ~ 10 10 years or during 
the cluster lifetime. 

One of the principal questions in understanding the 
long term evolution of cluster gas is whether mass as 
well as energy must be transported out from the central 
AGN. To investigate this possibility, we recently consid- 
ered idealized "circulation flows" in which gas heated 
near the center is buoyantly transported to large radii 
in the flow where it merges with the ambient gas and 
subsequently flows back in (Mathews et al. 2003; 2004). 
In these circulation flows gas moves simultaneously in 
both radial directions - there is no net flow of matter 
and no radiative cooling to low temperatures (satisfy- 
ing CR1). The dense inflowing gas emits most of the 
X-ray emission and resembles a traditional cooling flow 
that fills only a fraction of the available volume. The 
remaining volume is occupied by outflowing low-density 
bubbles of heated gas that do not contribute much to the 
X-ray spectrum. A single-temperature interpretation of 
the emission from such a flow is dominated by the inflow- 
ing gas and resembles a normal cooling flow (satisfying 
CR2). Furthermore, we showed that successful circu- 
lation flows must carry both mass and energy out from 
the center to approximately the cooling radius. The iron 
abundance peak that surrounds the central E galaxy is 
easily explained with circulation flows as the repository 
of all the iron produced by Type la supernovae in the 
central galaxy over time (satisfying CR3). 

Other heating scenarios can also be imagined. In the 
discussion below we explore the possibility that weak dis- 
sipating waves created near the center of the flow propa- 
gate throughout the hot gas, globally heating the gas and 
quenching the radiative cooling to low temperatures. 

The deep 200ks Chandra observation of the Perseus 
Cluster by Fabian et al. (2003) revealed the presence 
of an ensemble of approximately spherical weak shock 
waves ("ripples") within about 50 kpc of the center. The 
Perseus cluster also contains at least four X-ray cavities 
(bubbles) all within about 35 kpc of the center of the clus- 
ter (Birzan et al. 2004). Forman et al. (2003) describe 
similar spherical waves at 14, 17 and 37 kpc in the Virgo 
cluster. These distant shock-ripples are almost certainly 
created by the slow, intermittent inflation of new bub- 
bles near the central AGN. It is important to recognize 
that weak shocks are produced even by bubbles that ex- 
pand subsonically. The ripples observed have progressed 
about half way to the cooling radius which is about 130 
kpc in Perseus and 75 kpc in Virgo (Peres et al. 1998), 
assuming H n = 70 km s _1 Mpc -1 . 

Fabian et al. (2003) argue that dissipating shocks (rip- 
ples) excited by inner bubbles are likely to be the dom- 
inant heating mechanism that balances radiative cool- 
ing throughout the hot gas. They claim that alterna- 
tive modes of heating - large numbers of tiny efferves- 



cent bubbles (Begelman 2001), major central eruptions 
(Soker et al. 2001; Kaiser & Binney 2002) or heating 
near the buoyant bubbles (Churazov et al. 2001, 2002; 
Quilis et al. 2001; Briiggen et al. 2002; Bruggen & Kaiser 
2002; Bruggen 2003) - are less effective within r coo i than 
the quasi-continuous large scale dissipative heating by 
shock-ripples, which can persist even in those clusters 
(e.g. A1835) that do not currently have central bubble 
activity. 

The idea that ripple wave heating is dominant has been 
amplified recently by Clarke et al. (2005) who claim that 
"these shocks dissipate energy which may be sufficient 
to balance the effects of radiative cooling in the cluster 
cores." This conclusion has been strongly advocated by 
Ruszkowski, Bruggen & Begelman (2004a, b) and Dalla 
Vecchia et al (2004), who argue that the approximately 
spherical dissipation of large-scale waves is more efficient 
in globally balancing the radiative losses in cluster gas 
within r coo i than viscous heating associated with the flow 
just around the X-ray cavities. They also emphasize the 
importance of AGN intermittency in exciting individual 
ripple waves. 

Even with the best available observations, it is impos- 
sible to detect an increase in the gas temperature or ther- 
mal energy across individual ripples. The two ripples ob- 
served in Perseus, for example, are insufficient to increase 
the cluster energy (or entropy) enough to shut down the 
global cooling for many Gyrs. Because the amount of 
energy carried and entropy dissipated by each wave is so 
small, it is necessary to perform a calculation such as we 
describe here to test the idea that a large number of rip- 
ples sustained over time can provide enough dissipation 
to stop the cooling - without destroying the temperature 
profile as observed. We show that wave-heating fails be- 
cause the central regions of the cluster, through which all 
outgoing waves must pass, receive too much dissipated 
energy and after a few Gyrs become much hotter than 
the temperature observed. 

It is not surprising that Dalla Vecchia et al (2004) and 
Ruszkowski, Bruggen & Begelman (2004a, b) were able 
to reduce radiative cooling by ripple heating in their 2D 
and 3D simulations. To achieve this it is only necessary 
to adjust the (highly uncertain) power delivered to the 
hot gas by the AGN to be comparable to the bolomet- 
ric X-ray luminosity L x . By this ad hoc means radiative 
cooling can be shut down for a time ~ t coo i(rheat) within 
the radius r^ eat < r coo i where the AGN energy is de- 
posited by ripple shock waves. However, it is essential to 
determine whether or not cooling flows can be suppressed 
by ripple-heating over many Gyrs and still satisfy the gas 
temperature profiles generally observed (CR2). This is 
the question we address here. To establish with certainty 
that ripple wave dissipation can balance radiative losses 
in a cooling flow, it is necessary to calculate for times 
comparable to the relevant cooling time or evolutionary 
lifetime of the cluster. 

1.1. Viscous Heating by Shocks, Sound Waves and 
Turbulence 

There has been some confusion about the magnitude 
of the viscosity that is required for viscous dissipation 
to offset radiative losses in the hot gas in groups and 
clusters. In the absence of magnetic fields, the viscosity 
ft = 1.1 x iq- 16 T 5 / 2 g cm^ 1 s _1 depends only on the 
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temperature (Braginskii 1958). Viscous heating depends 
strongly on the local velocity gradient, ~ p,(du/d£) 2 erg 
cm~ 3 s~ 3 where i is a spatial coordinate appropriate 
to the scale of the region where dissipation occurs. Al- 
though viscous dissipation occurs on small spatial scales 
comparable to the mean free path, it can heat a large 
volume of gas by the passage of dissipating wave fronts 
or by widespread turbulence. 

X-ray cavities can heat the ambient gas by performing 
work and creating waves as they are initially formed and 
by viscous dissipation of turbulence produced by their 
buoyant motion. Birzan et al. (2004) estimated the 
PV work done by expanding cavities on the hot gas in 
16 galaxy clusters and concluded that "this mechanism 
alone probably does not provide a general solution to the 
cooling problem." If X-ray cavities inflate adiabatically 
into an invisid gas, an amount of work PV is transferred 
to the local gas. However, a significant fraction of this 
work is expended in displacing local gas outward in the 
cluster potential, but this gas returns to its original ra- 
dius after the buoyant bubble moves away, performing 
a nearly equal amount of negative work, — PV. Conse- 
quently PV overestimates the total energy received by 
the ambient gas from the cavities. 

The absence of resonance scattering of X-ray lines led 
Mathews et al. (2001) and Churazov et al. (2004) to 
suggest considerable subsonic turbulence in the central 
~ 60 kpc of Perseus but strong, extended turbulence 
may be inconsistent with steep iron abundance gradi- 
ents observed in Perseus and many similar clusters (e.g. 
Fabian et al. 2003). If such turbulent energy exists, the 
buoyant motions of X-ray cavities are a likely source. It 
is interesting to compare the rate of viscous dissipation 
near a cavity moving at ~ 30 kpc in Perseus (T « 3 x 10 7 
K; n e i=a 0.03 cm~ 3 ) to the local volume radiative emis- 
sivity. The drag force on a bubble (cavity) of radius Tb 
moving at velocity Ub is Fd = 5bpu 2 irrl where 5b < 1 is 
the drag coefficient. If the power expended against drag 
by a single bubble is FdUb, the total power within radius 
r in the flow is SbpufTrr 2 f 2 where f a — Nbirr 2 / Airr 2 < 1 
is the area filling factor for AT b bubbles at radius r in 
the flow. The total power deposited in volume 47rr 3 /3, 
(3/4)5bpulf 2 /r, is ultimately dissipated by viscosity at 
a rate 

ivisc = S b M 3 f a l^ = 2 ' 6 X 10 " 26 ^Va 2 



3 x 10 7 K 



3/2 



n e 



erg 



,0.03 cm- 3 / V30 kpc, 

where we assume that the bubbles move subsonically, 
i.e., M. — Ub/c s < 1 where c s = (■ykT/0.61m p ) 1 ^ 2 is the 
sound speed. Note that the absolute value of the viscosity 
determines the small scales of energy dissipation but not 
the amount of turbulent energy dissipated. The rate that 
energy is radiated at this radius is 

Srad= ( —) A(T,z)«2.4xlO- 26 



n e 



,0.03 cm" 



A 



2 x lO" 23 



erg 



Since we expect S b M 3 f 2 <C 1, the turbulent viscous dis- 
sipation directly associated with bubble buoyancy is un- 
likely to equal the local radiative losses. 



Viscous dissipation p(du/d£) 2 and heating can also 
be important in large amplitude sound waves and weak 
shock waves and must be investigated as a possible 
means of heating cluster gas. When ripple-waves were 
discovered in Perseus and Virgo, it was assumed that 
they were weak shocks of the sort that accompanies the 
creation of X-ray cavities (e.g. Brighenti & Mathews 
2002b). Recently Fabian et al. (2005) propose that clus- 
ter gas can be heated by linear sound waves of frequency 
/ ~ (3 x 10 5 yr)- 1 that dissipate on scales 



hi. 



2pc 3 



1 



(4 M /3) (27T/J2 



50 kpc. 



Fabian et al. show that it is possible to approximately 
offset radiative cooling in large clusters if the wave fre- 
quency is sufficiently high and the sound wave luminosity 
at this frequency is sufficiently large. The low amplitude 
linear waves required for this type of heating probably 
cannot be detected by X-ray observations. By compar- 
ison, in shock waves the dissipation occurs over a few 
mean free paths regardless of the magnitude of the vis- 
cosity, i.e. unlike linear wave dissipation the shock jump 
conditions and entropy increase do not explicitly depend 
on the viscosity. 

An important question that has not been throughly 
considered is the competition between the viscous dissi- 
pation of large amplitude sound waves as they decay with 
no change in the wave form and the progressive steepen- 
ing of the wave profile that results in a shock wave. The 
dimensionless ratio that determines these quite different 
evolutionary changes in the wave form is the Goldberg 
number 

(7 + l)M/2 



r = 



where Su = Mc s is the velocity amplitude of the wave 
and k = 2tt/X (e.g. Kinsler et al. 2000). The dimension- 
less Goldberg number is the ratio of M, an indicator of 
the degree of non-linearity in the wave, to t^ iss /k, the 
amplitude attenuation over one wavelength at the fun- 
damental frequency of the wave form. The evolution of 
waves created by X-ray cavities in the Perseus cluster 
with period T p can be estimated by combining the vari- 
ous expressions above, 



r = 350^ 



0.03cm- 



3 x 10 7 yr 



3 x 10 7 K 
f 



where = p/p s Ss 0.2 is the deviation of the viscos- 
ity from the Braginskii-Spitzer value. Since M. is less 
than (but comparable to) unity for subsonically expand- 
ing cavities, it is obvious that T ^> 1 so the waves created 
by expanding X-ray bubbles must steepen rapidly into 
weak shocks. 

1.2. Density Profile and Shock Wave Dissipation 

We describe below how hot cluster gas is heated by dis- 
sipating waves created near the central AGN. Since the 
gas density in groups and clusters varies less steeply than 
p oc r~ 2 , expanding ripple waves created in the bubble 
region (~ 10 — 20 kpc) lose amplitude as they encounter 
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an increasing mass per unit area; this loss in wave ampli- 
tude occurs in addition to viscous dissipation. We wish to 
determine if the temperature profiles typically observed 
in clusters are consistent with ripple-shock heating to ra- 
dius Theat where t coo i(rh e at) is comparable to the likely 
age of the cluster. Specifically we study wave dissipation 
in the Perseus cluster, where weak ripple-shocks are well- 
documented by observation. We begin with a discussion 
of the properties of this cluster. 

2. THE PERSEUS CLUSTER 

Churazov et al. (2004) give analytic fits to the density 
and temperature profiles in Perseus at the spatial res- 
olution of their XMM-Newton observations. The deep 
Chandra observation of of Perseus by Sanders et al. 
(2004) is consistent with the same temperature profile 
as Churazov et al., 

* [1 + (r/?1)3] K (1) 



8.12 x 10' 



[2.3 + (r/71)3] 

where r is in kpc with a distance determined with Hq = 
70 km s _1 Mpc -1 . The azimuthally averaged density 
profile in Perseus is given by 

0.0192 

n e (r) 



0.046 



[1 + (r/18) 3 ] 
0.0048 



[1 + (r/57) 2 ] 1 - 8 [1 + (r/200) 2 ] 



210.87 



cm 



(2) 



The last two terms are taken from Churazov et al. (2004) 
and the first term corrects for the Chandra observations 
of Sanders et al. in the high density core. 

The acceleration of gravity g(r) in the Perseus cluster 
can be found by assuming hydrostatic equilibrium in the 
hot gas and differentiating the equations above, 



9nT{r) 



kT ( 1 dn 



firrip 



l_eff\ 

n e dr T dr J 



(3) 



and this is shown with the dot-dashed line in Figure 1. 
Within about 10 kpc the observed acceleration g n T(r) 
decreases since T(r) and n e (r) in Equations (1) and (2) 
become constant as r — > 0. A centrally decreasing grav- 
itational acceleration is clearly inconsistent with the ex- 
pected dark and luminous mass distributions in this re- 
gion. The azimuthally (i.e. spherically) averaged fitting 
functions for n e (r) and T(r) fail for r ^ 10 kpc since this 
is the region occupied by large X-ray cavities. Even if the 
gas is close to hydrostatic equilibrium everywhere, as we 
believe to be the case, azimuthal averaging underesti- 
mates n e and d\ogn e /dr in the high density intercavity 
gas that produces most of the X-ray emission. Although 
the drop in g n T(f) within 10 kpc is unphysical, in the 
region 10 r 5s 300 kpc, where both n e (r) and T(r) 
are observed, g n r should accurately measure the Perseus 
gravitational field. 

We match the observed g n T(r) by combining an NFW 
dark halo and a stellar contribution from the cluster- 
centered galaxy NGC 1275, g = g n fw + 9*- The halo 
acceleration is simply 



9nfw — 



GM vir log(l + y) - y/(l + y) 



(4) 



r 2 log(l + c) - c/(l + c) 

where y — cr/r vir , c is the concentration, and r is the 
radius in kpc. The acceleration g n f w ( r ), determined with 



parameters M vir = 8.5 x 10 14 M Q , r vir = 2.440 Mpc, and 
c = 6.81, and shown as a long-dashed line in Figure 1 is 
a reasonably good fit to g n T{r) at 10 Ss r Ss 300 kpc. 

The physical properties of the central galaxy in 
Perseus, NGC 1275, are problematical since it appears 
to be an elliptical galaxy that has experienced a re- 
cent merger, perhaps with a smaller spiral galaxy (Con- 
selice et al. 2001). The global color of NGC 1275, 
(B - V) = 0.72 ± 0.03, is significantly bluer than a 
normal E galaxy, but the stellar surface brightness pro- 
file follows a de Vaucouleurs R 1 ^ 4 law to at least 150 kpc 
(Prestwich et al. 1997). We therefore assume here that 
the stellar mass of NGC 1275 is dominated by an old stel- 
lar population, as in normal E galaxies, intermixed with 
an additional population of young, luminous stars that 
does not contribute significantly to the total mass. The 
stellar mass distribution M*(r) and g»{r) for a de Vau- 
couleurs profile can be found from the effective radius, 
the total luminosity Lb of old stars, and the stellar mass 
to light ratio M/Lb- At a distance of 73.4 Mpc the ef- 
fective radius of NGC 1275 is R e = 6.41 kpc (RC3). The 
total magnitude My = —21.62 of the old stars can be 
found from the stellar velocity dispersion log er* = 2.40 
(Bettoni et al. 2003) using the log cr* — My correlation 
from Faber et al. (1997). If B — V = 0.95, as in normal E 
galaxies with old stellar populations, then Mb = —20.67 
and Lb = 2.70 x 10 10 L,b©, and we ignore the mass and 
luminosity of recently formed stars. For this Lb we find 
M/L B = 9.0 from Trujillo, Burkert & Bell (2004) so the 
total stellar mass of NGC 1275 is M H = 2 A3 x 10 11 M©. 
The stellar acceleration g* (r) = 67 M* (r) / r 2 for a de Vau- 
couleurs profile, shown as the short-dashed line in Figure 
1, can be accurately fit with 

-l/s 



9*{r) = 



„0.5975 



3.206 x 10" 



+ 



1.861 x 10" 6 



in cgs with s — 0.9 and r in kpc. We neglect the self 
gravity of the hot gas. The total two-component accel- 
eration g(r) = g n fw +9*, shown as a solid line in Figure 
1, fits gnrir) quite well for r > 10 kpc. 

3. EXCITING WAVES 

We wish to generate spherical shock waves of low am- 
plitude that propagate out through the Perseus gas. The 
objective is to determine if this type of shock heating 
(dissipation) can provide enough energy to the cluster 
gas to shut down the cooling for times <~ 4 Gyrs com- 
parable to the likely age of the cluster and still preserve 
the observed temperature profile. Waves can be gener- 
ated with a sinusoidally oscillating piston located in the 
bubble region. To create such a AGN wave machine it is 
convenient to use Lagrangian rather than the usual Eule- 
rian hydrodynamics since the boundary condition at the 
piston is more straightforward. The equation of motion 
is 

+ «-< 

where 

dr 



u = 



dt' 



(6) 



is the velocity defined in terms of the local Eulerian ra- 
dius r. The equation of continuity is simply 

pA(47rr 3 /3) = Am (7) 
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where A(47rr 3 /3) is the volume of computational zones 
of mass Am. The thermal energy equation with radiative 
losses is 



de 
dt 



dp 



dt 



_P_ 

rrir 



A(T,z) 



(8) 



where e = 3P/2p is the specific thermal energy and z w 
0.4 is the metal abundance in solar units. The artificial 
viscosity 

Q = a 2 p{Auf (9) 

depends on the velocity difference Am across the compu- 
tational zones and a dimensionlcss coefficient a of order 
unity to smooth the post-shock flow. The artificial vis- 
cosity ensures that the correct dissipative heating in the 
shock occurs over several computational zones. 

Spherical ripple waves can be generated by requiring 
that the inner boundary r p (t) of the innermost compu- 
tational zone, initially at r p0 , oscillates with amplitude 
Ar p and period T p , 



r p = r p0 + Ar p sin(27rf/T p ) 



(10) 



with velocity 



u p = (2vrAr p /r p )cos(27rt/r p ). 

Gas within the initial piston radius r p o does not partic- 
ipate in the flow - this is reasonable since we are inter- 
ested in heating the much larger volume of gas between 
the 10-30 kpc region where the cavities are observed and 
the distant cooling radius, r coo i w 130 kpc. The total me- 
chanical energy delivered to the gas by the piston after 
time t is 



E, 



,(*) = f ^r 2 p P p u p dt 
Jo 



(11) 



where P p (t) is the instantaneous gas pressure at r p (t) . If 
P p were constant then no mechanical luminosity (power) 
is delivered to the gas in one cycle, 



(L m )cy — (dE m /dt) 



cy 



8Tr 2 r 2 pa Ar p P p 



j [1 + (Ar p /r p0 ) sinr] 2 cosrdr = 0. 
Jo 



The mechanical luminosity depends critically on the 
asymmetry in the gas pressure at the piston P p (t) which 
is significantly larger when the piston moves outward, 
compressing the local gas, than when it recedes inward. 

Ar p /T p oc 



' p0 



Nevertheless in general we expect L Tl 
r 2 \u p \. The piston parameters - r p o, Ar p , and T p - can 
be selected by choosing T p ~ 10 7 — 10 s years, as ex- 
pected for the AGN periodicity from many observations, 
and then selecting Ar p so that the peak piston velocity 
| Up | = 2itAr p /T p remains subsonic since strong shocks 
are not observed around the X-ray cavities. Finally r p o 
can be increased until (L m ) is comparable to the bolo- 
metric luminosity L x , a necessary condition to quench 
the cooling flow. Alternatively, we can compare E m (t) 
with the total energy radiated after time t, 



E rad (t) = dt (p/m P YA(T, z)dV, 
Jo Jv 

where V is the total computational volume. 



4. A PURE COOLING FLOW IN PERSEUS 

Before considering the heating effects of outwardly 
propagating waves, it is useful to calculate a pure cool- 
ing flow in the Perseus cluster environment. This can be 
accomplished by fixing the piston at some small radius 
r p o = 1 kpc and setting Ar p = 0. In our solutions of the 
Lagrangian equations we begin with a configuration that 
is in strict hydrostatic equilibrium, which is probably an 
excellent approximation outside the bubble region and a 
good approximation even near the center where the gas 
velocities are largely subsonic. Initial hydrostatic equilib- 
rium is enforced on the computational grid by demanding 
that the differenced form of equation (5) corresponds to 
zero velocity Uj — at all zone boundaries rj, 
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9j 



k (p 



0.5(p j+ i /2 + Pj-i/ 2 ) 



j + l/2^ + 1/2 



Pj-l/2^-1/2) 



(12) 

pm p 0.5(rj+i - rj-i) 

The gravitational acceleration at every zone boundary is 
<7j = g(rj), where g(r) is plotted in Figure 1 with a solid 
line. The observed temperature profile is insensitive to 
the volume filling factor in the bubble region (Mathews, 
Brighenti & Buote 2004) and Chandra and XMM obser- 
vations are consistent with the same T(r), so it is natu- 
ral to assume Tj+1/2 = T{rj +1 / 2 ), using the analytical fit 
to the observed temperature profile described previously. 
With these assumptions, the difference equation above 
can be solved recursively for Pj +1/2 from the innermost 
zone outward. If necessary, the density of the innermost 
zone is adjusted slightly so that the overall density profile 
agrees with that described by equation (2). In general 
the initial gas density Pj+1/2 within r ~ 10 kpc, which 
corresponds to the intercavity density, is somewhat larger 
than the analytic fit (Eqn. 2), which is a spherical aver- 
age including both cavity and intercavity gas. 

As gas radiates and flows toward the center in a nor- 
mal cooling flow, the density rises until radiative cooling 
eventually becomes catastrophic. When the gas temper- 
ature in a Lagrangian zone drops below 3 x 10 5 K, the 
zone becomes spatially very narrow and we remove it 
from the calculation, filling its small volume with a small 
amount of gas from neighboring zones. Since the inner 
boundary at r p o = 1 kpc is stationary, all the gas cools 
just beyond that small radius. A shock near r p o prevents 
the inflowing gas from attaining large negative velocities. 

The solid lines in Figure 2 show the cooling flow after 
4 Gyrs compared to the initial profiles shown as dashed 
lines. We select 4 Gyrs (redshift z ~ 0.5) since massive 
clusters such as Perseus may have evolved significantly 
over longer times. Toward the center of the flow the 
computed gas density n e rises above the observations in 
r < 20 kpc as r — > r p0 ; such a density rise is typical in 
traditional cooling flows (see e.g. Mathews & Brighenti 
2003) and is enhanced further in Perseus by the central 
X-ray cavities that lower the mean observed density. By 
4 Gyrs the gas cooling rate shown in Figure 3 has sta- 
bilized at M = 250 M Q yr -1 and the huge total mass 
of cooled gas at the center, M coo ; = 9.8 x 10 11 M Q , is 
equal to all the gas within 73 kpc in the initial configu- 
ration. Our M is nearly identical to the cooling flow rate 
for Perseus predicted by Sanders et al. (2004), M w 255 
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M yr _1 , but is inconsistent with their Chandra obser- 
vations that show no evidence that gas is cooling below 
~2x 10 7 K. Heating by ripple wave dissipation or other 
means is necessary to prevent this cooling. 

We have advocated that gas near the centers of cooling 
flows is heated and buoyantly recirculated to large radii 
extending to ~ r coo ; (Mathews et al. 2003; Mathews, 
Brighcnti & Buotc 2004). In these "circulation flows" 
both mass and energy are transported outward, but the 
thermal profile of of the flow is nearly identical to that 
of a normal cooling flow, i.e. T(r) is insensitive to the 
radial variation of the volume filling factor occupied by 
the rising cavities (bubbles) as shown by Mathews et 
al. (2004). For this reason it is noteworthy in Figure 2 
that the computed cooling flow temperature profile af- 
ter 4 Gyrs is almost identical to the observed profile for 
r > 30 kpc. Since a pure cooling inflow matches the 
observed temperature and density in this outer region 
of Perseus so well, we infer that there is no need for 
conductive heat transport in this region of the flow, al- 
though some modest heating may occur. In the following 
section we explore whether this radiative cooling in the 
Perseus cluster, M = 250 M© yr _1 , can be significantly 
reduced by dissipative losses in outwardly propagating 
weak shock waves (ripples). 

5. RIPPLE-HEATED FLOWS 

To create a series of outward propagating weak shock 
waves that move out into the Perseus atmosphere, we 
allow the spherical piston at radius r p o to undergo si- 
nusoidal oscillations with amplitude Ar p and period T p . 
Since strong shocks are rarely observed surrounding X- 
ray cavities, we select parameters for which the maxi- 
mum velocity of the piston |it p>rna;c | = 2irAr p /T p is less 
than the sound speed in the local hot gas, c s (r p ) — 
[5fcT(rp)/3^m p ] 1 / 2 = 958[T(r p )/3 x I0 7 K] 1 / 2 km s" 1 . 
Note that weak shocks are created even when |u p | < c s . 

The flows with piston parameters listed in Table I 
are a representative sample of our exploration of heat- 
ing by dissipating ripple waves. The first flow of this 
type is shown in Figure 4 at two times, designated al 
after 4 x 10 7 yrs and a after 4 Gyrs. If the gas density 
varied as p oc r~ 2 , outwardly propagating dissipation- 
less waves would encounter nearly the same mass per 
unit area at every radius and their amplitude would re- 
main unchanged. However, in Perseus the density drops 
slower than p oc r -1-6 so the wave amplitude is expected 
to decrease with radius as a larger mass participates in 
the waves. This is seen most clearly in the gas velocity 
panel in Figure 4 for model a. The piston parameters of 
model a have been selected so that the mechanical en- 
ergy E m (t) expended by time t is initially slightly larger 
than the radiated energy, as illustrated in Figure 5 - this 
requires a piston velocity \u p \ that is a large fraction of 
the initial sound speed. As can be seen in Figure 3, the 
total mass of cooled gas in model a is small and this 
cooling occurs close to the piston. Nevertheless, this so- 
lution is clearly unacceptable because it has the wrong 
density and temperature profiles. As expected, most of 
the wave dissipation occurs at small radii in the flow, 
resulting in temperatures exceeding 10 8 K with negative 
dT/dr continuing to nearly the cooling radius. The view 
of model a at an earlier time in the first column of Fig- 
ure 4 (model al), shows that rather strong shocks were 



present at early times. By 4 Gyrs these shocks produced 
the heating shown in the second column (model a). Af- 
ter the central regions in Perseus have been heated to 
these unrealistically high temperatures, the piston pro- 
duces a more rounded innermost wave profile in flow a, 
a subsonic non- linear wave that shocks further out. 

In model b we consider a flow with the same velocity 
| Up | as in model a but with a shorter period T p — 10 7 yrs 
and the resulting flow has an entirely different character. 
The outer part of flow b beyond 50 kpc in Figure 4 is 
essentially a cooling inflow, and this is verified in the 
bottom panel where (u) < 0. The cooling rate and total 
cooled gas in model b are unacceptably large (Figure 3). 
All of this cooling occurs at r « 50 kpc where the gas 
temperature has a sharp minimum and the density has 
a corresponding rise. The waves generated by the piston 
arc unable to penetrate the impedance barrier at r w 50 
kpc and this accounts for the very low wave amplitude at 
larger r. Within r « 50 kpc the dissipating waves reflect 
back and forth raising the gas temperature far above that 
observed in Perseus. 

In models c — / the piston amplitude Ar p and period 
T p are chosen so that the piston velocity \u p \ is more 
subsonic than the previous solutions (Table 1). As we 
have discussed already, to achieve an approximate bal- 
ance between mechanical and radiative luminosity when 
| Up | is smaller, it it necessary to increase the mean piston 
position r p0 . In model c the piston position r p0 = 31.5 
kpc is near the outer limits of the observed X-ray cav- 
ities in Perseus. In Figure 5 it is seen that E m is only 
slightly less than E rac [ during the entire calculation to 
4 Gyrs, but Figure 3 shows that very little gas (all near 
the piston) has cooled. However, as with model a the gas 
temperature in the region 37 r < 57 kpc of model c 
exceeds 10 8 K, which is far in excess of the temperature 
observed in Perseus. 

Finally, we consider a series of models, d — f in which 
| u p | remains fixed at the small value in model c, but 
r p o is slowly reduced. A smaller mean piston radius r p o 
is expected to result in lower mechanical energy input 
E m oc r 2 |u p | relative to the total radiated energy E ra d- 
This increasing discrepancy can be seen by comparing 
models c and / in Figure 5. The overall character of solu- 
tions c — ► e remains about the same, with a small amount 
of cooling (Figure 3) but with temperature profiles that 
peak at small radii and therefore are inconsistent with 
the observations. Finally, in model / the mechanical en- 
ergy is unable to keep large amounts of gas from cooling 
(Figure 3). Unlike the flows c — e, most of the cooling in 
flow / occurs not at the piston but in the cooling cusp 
at r w 30 kpc. Model / resembles model b in having 
an external cooling flow in r > 30 kpc that is supported 
by gas heated by waves that reflect at this impedance 
barrier. 

The impedance barriers in models b and f, where the 
density and temperature vary significantly over one wave- 
length, are a long term result of unrealistically intense 
central wave-heating, and are grossly incompatible with 
observations. However, long before impedance barriers 
appear in the solutions as shown in Figure 4, the tem- 
perature and density structure of the hot gas has evolved 
far from typically observed profiles. In principle the hot 
central cores in these models could be smoothed by out- 
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ward thermal conduction, but this is irrelevant since in 
observed clusters the inner temperature gradient is pos- 
itive and the conduction, if it were important, should be 
heating the central gas, not cooling it. The sharp posi- 
tive density jump in models b and f are very likely to be 
disrupted by Rayleigh- Taylor instabilities. However, in 
our previous 2D calculations, in which the RT instability 
was allowed to develop naturally, we show that includ- 
ing more computational dimensions does not save the 
problems associated with the wave-heating hypothesis 
(Brighenti & Mathews 2003). All centrally heated flows 
we considered (allowing for RT) resulted in dT/dr < 
in violation of CR2. 

6. SUMMARY AND CONCLUSIONS 

Following the discovery of shock-ripples in the Perseus 
and Virgo clusters, it has often been assumed that this 
is the fundamental global heating mechanism that pre- 
vents the gas from cooling as required by X-ray spectral 
observations. A heating process can dynamically shut 
down a cooling flow if it distributes enough energy to 
reduce or reverse the small negative flow velocity out to 
approximately the cooling radius. It is clear that the 
subsonic formation of large X-ray cavities creates weak 
shocks that move far out into the flow and that these 
waves dissipate some of their energy in remote regions of 
the flow. 

Nevertheless, we have demonstrated here that the 
transport of energy by outwardly propagating waves is 
not likely to be the dominant process that accounts for 
the absence of cool and cooling gas in these flows. The 
main problem with this idea is that the density struc- 
ture of the hot gas varies as p oc r~ s where s < 2 so 
that outgoing waves dissipate most of their energy in the 
central part of the flow and too little at r ~ r coo ;. The 
result is that after a few Gyrs the gas temperature flat- 
tens and eventually peaks near the center of the flow, just 
opposite to observed temperature profiles in cooling core 
clusters (e.g., Allen, Schmidt, & Fabian 2001), violating 
the second cardinal requirement (CR.2). Much enthu- 
siasm for ripple heating has resulted from calculations 
that proceeded for only ~ 0.3 Gyrs (e.g. Ruszkowski, 
M., Briiggen, & Begelman 2004a, b), as in our model al, 
but we have shown here that the temperature profile de- 
teriorates on more relevant Gyr timescales. In a similar 
calculation dalla Vecchia et al. (2004) describe wave- 
heated flows and compute to 1.5 - 5 Gyrs, but find, as 
we do, that the resulting temperature profile differs from 
those observed. 

Very recently, after this paper was submitted, Fujita & 
Suzuki (2005) presented a model of steady state cluster 
gas heating with weak shocks that reproduces many of 
the same results derived here. In particular they find that 
a series of weak shocks propagating out in the Perseus 
Cluster results in unrealistically strong heating in the 
core. 

The unsatisfactory results described here are similar 
to the difficulties we encountered in our earlier explo- 
rations of cooling flows that were heated in a variety of 
ways (Brighenti & Mathews 2002a; 2003). Our centrally 
heated cooling flows produced rising bubbles and weak 
shocks that propagated to large radii, but we did not 
discuss distant wave dissipation in detail. We also found 
that (feedback) heating the central regions of cooling 



flows often simply moved the cooling to larger radii but 
did not diminish the cooling rate (as in the ripple-heated 
models b and / described here). Large scale cooling in- 
flow can be stopped if the AGN feedback heating is rapid, 
extends to ~ r coo i and is supplemented by thermal con- 
duction from large radii (Ruszkowski & Begelman 2002), 
but the amount that thermal conduction is suppressed 
by magnetic fields needs some fine-tuning (Brighenti & 
Mathews 2003). A similar unphysical fine-tuning is ap- 
parent in the wave-heated models described here - the 
location, amplitude and period of wave creation by the 
X-ray cavities must be finely and unrealistically regu- 
lated to ensure just the proper amount of heating in dis- 
tant regions of the flow to balance radiative losses, i.e. 
E m ~ E ra< i, although the resulting temperature gradient 
is unacceptable. 

Most of this fine-tuning can be avoided if both mass 
and energy are redistributed outward in the hot gas as 
in the "circulation flows" described by Mathews et al. 
(2003) and Mathews, Brighenti & Buote (2004). In these 
flows the mass transported by dense inflowing gas is bal- 
anced by a nearly equal amount of gas flowing outward 
to large radii carrying both mass and energy. The denser 
inflow emits most of the radiation and appears as a nor- 
mal cooling flow with a temperature profile that matches 
the observations. Meanwhile, bubbles of hot gas heated 
by the central AGN flow out rapidly and contribute much 
less to the overall X-ray emission. Our previous 2D com- 
putational models of heated flows lacked the spatial reso- 
lution to isolate individually rising bubbles as they move 
large distances upstream through the denser gas, but cir- 
culation flows can be computed in a more schematic fash- 
ion (Mathews, Brighenti & Buote 2004). In mass circu- 
lating flows it is only necessary that gas near the central 
AGN be heated sufficiently to have an entropy similar 
to gas at approximately r coo i or beyond. The best argu- 
ment for outward mass circulation are the vast regions 
of iron-enriched gas that extend far beyond the stars in 
the central galaxy where the iron is produced. 

We argue here that neither the PV work done as the 
hot X-ray cavities inflate or viscous heating by cavity- 
driven turbulence is a strong source of local heating; this 
is consistent with the low gas temperatures observed in 
cooling-core clusters in the central region where bubbles 
are observed. This suggests that the AGN energy (and 
mass) must be distributed non-locally to more distant 
gas. 

To summarize, we have shown that viscous dissipation 
in outward propagating weak shock waves cannot keep 
the gas in the Perseus cluster from cooling without heat- 
ing the gas in the central regions far above the observed 
temperatures. Since weak shocks have been observed in 
Perseus and elsewhere, some heating is expected but they 
are not the dominant heating mechanism. The shocks 
observed are almost certainly produced by the subsonic 
inflation of X-ray cavities, which, as we have shown, pro- 
duce non-linear waves that rapidly steepen into weak 
shocks. Evidently cavity-produced shocks are too infre- 
quent or weak to generate mechanical luminosities com- 
parable to the X-ray losses. 

Studies of the evolution of hot gas in elliptical galax- 
ies at UC Santa Cruz are supported by NASA grants 
NAG 5-8409 & ATP02-0122-0079 and NSF grants AST- 
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TABLE 1 

PISTON PARAMETERS TO CREATE SHOCK-RIPPLES 



model 


t a 


r p 


Ar p 




M 


rl\u p \ 




(Gyr) 


(kpc) 


(kpc) 


(yr) 


(km/s) 


(kpc 2 km/s) 


al 


0.04 


23.8 


3.9 


3 x 10 7 


800 


4.5 X 10 s 


a 


4 


23.8 


3.9 


3 x 10 7 


800 


4.5 X 10 s 


b 


4 


23.8 


1.3 


1 x 10 7 


800 


4.5 x 10 s 


c 


4 


31.5 


2.23 


3 x 10 7 


457 


4.5 x 10 s 


d 


4 


23.8 


2.23 


3 x 10 7 


457 


2.6 x 10 s 


e 


4 


18.1 


2.23 


3 x 10 7 


457 


1.5 x 10 5 


f 


4 


13.8 


2.23 


3 x 10 7 


457 


0.87 x 10 5 



a Total time of flow calculation. 
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Fig. 1. — Acceleration of gravity in the Perseus cluster: g n T derived from the observed gas density and temperature (dot-dashed line), 
g n f w component from the dark halo (long dashed line), g t component from stars in the central galaxy NGC 1275 (short dashed line) and 
9 = 9nfw + 9* is the total adopted acceleration (solid line). 
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Fig. 2. — A standard cooling flow in the Perseus cluster after 4 Gyrs (solid lines) overplotted with the observed gas density, temperature 
and pressure (dashed lines) and initial velocity (dashed line). 
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Fig. 3. — Total mass of gas in Mq that cooled after time t in the cooling flow (model cf) and in six flows with ripple-shock heating, 
models a — > /. 
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Fig. 4. — Computed gas density, temperature, pressure and velocity after 4 Gyrs in models a — c and / (solid lines) compared with the 
observed (or initial) profiles (dashed lines). 
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Fig. 5. — Total amount of mechanical energy E m (t) delivered to the Perseus gas by the spherical piston after time t (dashed lines) for 
flow models a — c and / is compared to the total bolometric energy lost in radiation E ra d(t) (solid lines), both in ergs. 



